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Introduction 


In  the  last  few  years,  advances  in  computing  power  and  algorithm  development  have 
allowed  the  creation  of  highly  accurate  nonlinear  aero-elastic  simulation  codes  [1,2]. 
These  codes  typically  have  fluid  meshes  with  a  very  large  number  of  nodes,  and  a  lesser 
(but  still  large)  number  of  structural  nodes.  System  orders  in  the  millions  of  degrees  of 
freedom  are  not  uncommon.  Designing  control  systems  using  these  models  is  not 
feasible  with  current  techniques  due  to  these  large  sizes.  Instead,  some  method  or 
combination  of  methods  for  reduced  order  modeling  must  be  employed  to  create  a  low 
order  model  of  the  aeroelastic  system. 

The  most  obvious  model  reduction  method  to  try  on  a  large-scale  system  is  eigen 
truncation.  Reduced  order  models  based  on  eigenreduction  techniques  are  attractive  to 
the  control  designer  because  of  their  common  usage  in  structural  dynamics  control  (see 
review  in  [3]).  Many  nonlinear  solvers  have  built  in  parallelization  and  linear  equation 
solvers  than  can  be  adapted  to  solve  for  the  first  few  eigenvectors,  which  can  then  be 
used  as  a  basis  for  a  reduced  order  model.  The  current  paradigm  is  that  modal  truncation 
is  not  a  good  technique  to  use  in  aeroelasticity  because  the  fluid  eigenvalues  tend  to  be 
closely  spaced,  necessitating  a  large  number  of  retained  eigenmodes  to  achieve  accurate 
results.  Instead,  emphasis  has  been  placed  on  methods  that  map  system  inputs  to  system 
outputs  directly,  including  approximate  balanced  truncation  methods  [4]  and  the  Proper 
Orthogonal  Decomposition  (POD)  method  [5,6].  The  impetus  for  this  paper  grew  out  of  a 
desire  to  use  a  reduced  order  model  of  the  fluid  coupled  to  a  reduced  order  modal 
truncation  model  of  the  structure. 

Once  a  reduced  order  model  is  created,  a  controller  can  then  be  designed  about  it 
using  conventional  techniques.  Problems  can  arise  once  the  controller  based  on  the 
reduced  order  model  is  applied  to  the  full  order  system.  Since  the  controller  has  no 
knowledge  built  into  it  about  the  unmodelled  states  in  the  full  system,  interactions  can 
excite  those  states  into  instability.  Previous  work  using  modal  truncation  based  ROMs 
[7]  as  well  as  coupled  systems  of  modal  based  ROMs  [8]  has  shown  that  by  using  a 
Residual  Mode  Filter  (RMF)  to  filter  out  the  unwanted  interaction,  stability  can  be 
restored  to  the  system.  These  previous  stability  proofs  rely  on  the  modal  nature  of  the 
systems  and  are  not  valid  once  the  coupling  terms  created  by  balanced  truncation 
methods  are  considered. 


1  Residual  Mode  Filters  and  Actuated  Systems 


The  Residual  Mode  Filter  (RMF)  is  a  technique  that  can  be  used  to  regain  stability  in  a 
high  order  system  driven  by  a  reduced  order  controller  when  unmodelled  interactions 
have  caused  instability.  In  general,  RMF’s  work  by  splitting  the  state  vector  x  into 
components  xn,  xr,  and  xq  where  xn  contains  those  states  in  the  reduced  order  model 
(including  all  initially  unstable  states),  xr  contains  the  remaining  stable  modes,  and  Xq 
contains  those  modes  that  have  been  driven  unstable  through  feedback  control.  The  RMF 
is  then  applied  to  cancel  the  xq  modes  from  interacting  with  the  controller,  thereby 
causing  the  system  to  regain  stability. 
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Now  if  instead  of  a  single  dynamical  system  there  is  a  coupled  problem,  things 
become  more  complicated.  Utilizing  the  method  described  in  [9],  the  system  can  be  split 
into  a  coupled  aeroelastic  response  and  an  aerodynamic  actuation  dynamic.  In  this 
formulation,  one  system  can  be  seen  acting  as  an  actuator  on  the  other,  and  the  combined 
system  can  be  denoted  as: 

x  =  Ax  +  Bu 

y  =  Cx  +  Du\  jc(0)  =  xn 

w  0  (1) 

z  =  Fz  +  Huc 
u-Lz 

where  z  denotes  actuator  states  (aerodynamic  actuation  dynamic)  and  x  system  states 
(coupled  aeroelastic  reponse).  The  system  for  the  case  we  are  considering  can  then  be 
decomposed  into  modal  form: 

x^  A^Xfl  -f-  Bnu 

*r  ~  ArXr  +  Bru 

<  #  (2) 

Xq  —  AqXq  +  BqU 

y  ~  C N  xN  +  CR  xR  4-  Cq  Xq  +  Du 

where  Xn,  Xr,  and  xq  are  defined  as  before.  Note  that  all  the  inherently  unstable  modes 
are  collected  in  Xn,  therefore  Ar  and  Aq  will  be  stable  matrices  (in  the  sense  that  all  the 
eigenvalues  will  lie  in  the  left  half  plane).  An  output  feedback  state  estimator  controller 
designed  for  this  system  will  have  the  following  form: 

xN  =  Anxn  +  BnLz  +  Kn  {y  -  yN ) 

yN  =Cnxn+DLz 

i  =  Fz  +  Hue 

Uc  = 

Because  of  how  we  have  defined  the  states,  this  controller  applied  to  the  full  order 
system  will  be  stable  if  the  xq  terms  are  deleted.  If  the  following  error  terms  are  then 
defined: 


This  system  is  by  definition  stable.  Note  in  particular  the  stable  sub-block  consisting  of 
the  first  4  rows  and  columns  of  eqn.  (5),  we  will  be  using  this  sub-block  later.  If  we  next 
introduce  the  unstable  modes  xq  back  into  our  system,  some  method  to  compensate  for 
their  instability  is  required.  A  RMF  can  be  used  to  filter  the  output  signal  y  being  fed  into 
the  state  estimator  in  this  manner: 
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\xN=ANxN+BNLz  +  KN(y-yN-yQ) 
9n  =Cnxn+DLz 
z  =  Fz  +  Hur 


AqXq+BqU 

9q~  Cq*q 

Mc  ~  GnXn 


(6) 


By  introducing  a  new  error  term  eQ  defined  as: 


Xq  Xq 


(7) 


a  new  matrix  equation  can  be  written  that  now  incorporates  the  entire  system: 


'  Ah 

0 

0 

bnl 

0 

0 

0 

0 

(an-kncn) 

kncr 

0 

0 

-KnCq 

(bnl-kndl) 

0 

0 

brl 

0 

0 

0 

hgn 

hgn 

0 

F 

0 

0 

0 

0 

0 

0 

bql 

A? 

0 

0 

0 

0 

0 

0 

0 

aq 

bql 

0 

0 

0 

0 

0 

0 

F 

(8) 


Note  the  sub-block  formed  by  the  first  four  rows  and  columns;  this  is  the  same  sub-block 
as  in  eqn.  (5),  and  was  by  definition  stable.  Because  of  this,  the  stability  of  the  entire 
matrix  is  seen  by  successively  including  an  additional  row  and  column  of  eqn.  (8)  and 
noting  the  diagonal  nature  of  subsequent  subsystems.  Since  all  the  inherently  unstable 
modes  were  collected  into  An,  Aq  is  stable,  and  F  is  also  stable.  Therefore  the  entire 
matrix  given  in  eqn.  (8)  is  stable. 


2  Residual  Mode  Filters  with  Reduced  Order  Balanced  Actuators 

For  the  actual  problems  of  interest  in  this  work,  which  consist  of  fluid  systems  acting  as 
actuators  and  coupled  to  structural  systems,  it  is  not  practical  to  include  the  entire  model 
of  the  fluid  system  within  the  controller.  If  a  ROM  of  the  fluid  is  used  that  is  based  on  a 
balanced  reduction  method,  the  actuator  can  be  represented  by: 

^■N  =  FnZ-N  FnrZr  +  H Nuc 

‘  zR  =  FRzR  +  Frnzn  +  H Ruc  (9) 

u  —  Lnzn  +  Lrzr 

where  zn  represents  those  states  that  are  retained  in  the  ROM  and  zR  those  states  that  are 
discarded.  Note  that  the  key  difference  to  the  work  described  in  [8]  is  the  inclusion  of  the 
cross-coupling  terms  Fnr  and  Frn.  Because  of  the  nature  of  a  balancing  reduction,  the 
influence  the  zR  states  have  on  the  output  u  and  the  retained  states  zn  is  small  (and  can  be 
made  smaller  by  including  more  states  in  zn).  Because  of  this,  it  is  natural  to  introduce  a 
singular  perturbation  term  in  front  of  the  zR  equation,  and  analyze  what  happens  as  this 
term  is  varied  from  zero  to  some  small  positive  value: 
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(10) 


Z-N  ~  BNZN  ■*"  FnrZr  +  H NUc 
•  £ZR  =  Frzr  +  Frnzn  +  Hruc 
u  —  Lnzn  +  Lrzr 

When  e  =  0,  zr  can  be  solved  for  in  terms  of  zn  and  uc: 

Zr  =~Fr  fr„Zn -F~'Hruc 

zn=Fnzn+Fnr(-F~'Frnzn -F~'Hruc)+Hnuc  (11) 
u  =  Ln zn  +  Fr (—  Fr  Frn zn  —  F *  H Ruc) 

If  the  following  substitutions  are  then  made: 


fr=-lrf;'hrgn 

F'  =  -/  F~'F 

1  R  ^R1  R  1  RN 


(12) 


Then  the  actuation  equations  can  be  written  as: 

Zn  =  {^n  ~  FNR Fr  Frn  )zn  +  {h  n  ~  FNRFR  H R  )Gwxw 
'  £Zr  ~  Frzr  +  Frnzn  +  H Ruc 
u=(ln  +  Fr)zn  +  Frxn 

Based  on  this  definition  of  the  actuator,  the  control  equations  become: 
\xN  —  An xn  +  Bn Ln zn  +  Kn  [y  —  yN  ~  9 q) 

A  /—I  A 

yN  ~~  fnxn 

Zn  =  i^N  —  Fnr  Fr  Frn  )zw  +  [h N  —  FnrFr  H r  )gnxn 
xq  ~  +  Fq  (ln  +  Fr  )zw  +  BqFrxn 


(13) 


(14) 


yQ  =  cnx 


Q^Q 


Uc  =  GNXN 


The  system  (defined  to  be  stable)  with  the  xq  modes  deleted  and  e  =  0  is  given  by: 

(15) 


d>=  Anco 


Where: 


0)  = 


Zn 

e. 


(16) 


And: 


A, 


( An  +  Bn  Fr  ) 

bnfr 

0 

bn{ln  +  f') 

0 

(~bnfr) 

4 

FnCr 

~BnK 

(BnLn) 

BrFr 

BrFr 

Br{fn+F') 

0 

z, 

Z, 

0 

z2 

0 

0 

0 

0 

0 

z2 

(17) 


With: 
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e,=an-kncn-bnfr 

■z,=(«„- os) 

^2  ~  Fn  ~  Fnr  Fr  Frn 

Proving  that  the  system  is  stable  with  the  Xq  modes  added  back  in  and  an  epsilon 
not  equal  to  zero  requires  the  use  of  the  Klimushchev-Krasovskii  theorem  [10], 
formulated  for  the  following  equation: 

(d)  —  Auco+  Ai2zR 
=  A2X(D+  A22zr 

Given  a  linear  singularly  perturbed  system  of  the  form  in  eqn.  (19)  with  A22  andAu- 
A12A22 1  Az\  stable,  then  there  exists  some  £o>0  such  that  for  each  0<£<£o,  eqn.  (19)  is 
stable. 

We  will  use  this  lemma  by  writing  the  complete  system  in  the  form  of  eqn.  (19), 

with: 


The  matrices  An,  A\2,  A21,  and  A22  are  then: 


(.Aw  +  BnFr) 

BnFr 

0 

bn{ln  +  f;) 

0 

0 

0 

(- bJr ) 

£* 

KnCr 

-BnF' 

0 

KnCq 

{bJr) 

BrFr 

BrFr 

Ar 

Br{ln  +  f;) 

0 

0 

0 

z, 

Z, 

0 

Z2 

0 

0 

0 

bqfk 

BqFr 

0 

bq{ln  +  f') 

Aq 

0 

0 

0 

0 

0 

0 

0 

Aq 

Bq  Fn 

0 

0 

0 

0 

0 

0 

z2  J 

\"„Gn 


000] 


A'n  "  Fr 


To  use  the  K-K  theorem,  we  need  to  show  that  Ai  1  is  stable  and  that  Ah-  A12A22" 

*  A2i  is  stable.  The  stability  of  An  can  be  seen  from  examining  An  in  eqn.  (17)  and  then 
applying  the  successive  partitions  as  shown  in  eqn.  (21),  since  Aq  is  defined  to  be  stable 
and  Z2  is  shown  to  be  stable  from  eqn.  (17).  An  -  Ai2A22]A2i  is  also  stable  since  it  equals 
Ai  1 .  Therefore,  the  entire  system  is  stable  due  to  the  K-K  theorem  for  some  epsilon  0  <  £ 

<e0- 
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3  Results 


This  proof  followed  the  outline  for  a  modal  based  actuator  given  in  [8],  where  an 
example  using  an  Aeroelastic  Piston  model  was  given,  which  is  based  on  the  ideas 
developed  in  [1 1].  The  same  model  was  used,  but  now  with  a  reduced  order  balanced 
model  of  the  actuator.  An  output  feedback  state  estimator  controller  was  designed  using 
the  reduced  order  system  and  actuator  that  gave  good  results  when  applied  to  the  ROM. 
When  the  same  controller  was  applied  to  the  FOM,  the  interaction  with  an  unmodeled 
mode  drove  the  response  unstable.  Including  a  RMF  as  in  eqn.  (14)  causes  the  full  order 
response  to  be  stabilized,  as  seen  in  Figure  1 . 


A.  Piston  Velocity  -  Reduced  System  with  Output  Feedback  Control 


V  v^'"' 

.V 

- 

0  0.05  0.1  0,15  0.2  025  0.3  0.35  04  0.45  0.5 

Time  (sec.) 


Figure  1 :  Residual  Mode  Filter  Example 


References 

[1]  P.  S.  Chen,  P.  S.,  Farhat,  C.,  Lesoinne,  M.,  and  Lanteri,  S„  Parallel  Heterogeneous  Algorithms  for  the 
Solution  of  Three-Dimensional  Transient  Coupled  Aeroelastic  Problems,  AIAA  36th  Structural 
Dynamics  Meeting,  New  Orleans,  LA,  1995. 

[2]  Pak,  C.,  Friedmann,  P.,  &  Livne,  E„  Transonic  Adaptive  Flutter  Suppression  Using  Approximate 
Unsteady  Time  Domain  Aerodynamics,  32nd  AIAA/ASME/ASCE/AHS/ASC  Structures,  Structural 
Dynamics  and  Materials  Conference ,  Baltimore,  MD,  1991. 

[3]  M.  Balas,  Trends  in  Large  Space  Structure  Control  Theory:  Fondest  Hopes,  Wildest  Dreams,  IEEE 
Trans.  Auto.  Controls,  AC-27(3),  June  1982. 

[4]  Baker,  M,,  Mingori,  D.L.  and  Goggin,  P.J.,  Approximate  Subspace  Iteration  for  Constructing 
Internally  Balanced  Reduced  Order  Models  of  Unsteady  Aerodynamical  Systems,  AIAA  Paper  #  96- 
1441, 37th  AIAA/ASME/ASCE/AHS/ASC  Structures,  Structural  Dynamics,  and  Materials  Conference, 
Salt  Lake  City,  UT,  Apr.  1 5- 1 7,  1 996. 

[5]  Wilcox,  K.  and  Peraire,  J.,  Balanced  Model  Reduction  Via  the  Proper  Orthogonal  Decomposition, 
AIAA  Paper  #  2001  -26 1 1 ,  1 5th  AIAA  Computational  Fluid  Dynamics  Conference,  Anaheim,  CA,  1 1  - 
14  June,  2001. 


8 


[6]  Mortara,  S.A.,  Slater,  J.C.,  and  Beran,  P.S.,  A  Proper  Orthogonal  Decomposition  Technique  for  the 
Computation  of  Nonlinear  Panel  Response,  AIAA  Paper  #  2000-141 1,  AIAA  Structures,  Structural 
Dynamics  and  Materials  Conference ,  Atlanta,  GA,  2000. 

[7]  Balas,  M.  Finite  Dimensional  Controllers  for  Linear  Distributed  Parameter  Systems:  Exponential 
Stability  Using  Residual  Mode  Filters,  Journal  of  Mathematical  Analysis  and  Applications ,  133(2), 
Aug.  1988. 

[8]  Hindman,  C„  Balas,  M.  and  Lesoinne,  M.,  Exponential  Stability  of  Controllers  Using  Reduced  Order 
Actuator  and  System  Models  with  Residual  Mode  Filters,  2000  American  Control  Conference , 
Chicago,  EL,  2000. 

[9]  Lesoinne,  M.  and  Hindman,  C.,  A  Mixed  Model  Reduction  Technique  for  Aeroservoelasticity,  AIAA 
Paper  #  99-1266,  40th  A1AAJASCE/AHS/ASC  Structures,  Structural  Dynamics,  and  Materials 
Conference,  St.  Louis,  MO.,  Apr.  12-15,  1999. 

[10]  Vidyasagar,  M.,  Nonlinear  Systems  Analysis,  Prentice-Hall,  N.J.,  1978. 

[1  l]Piperno,  S.,  Larrouturou,  B.  and  Lesoinne,  M.,  Analysis  and  Compensation  of  Numerical  Damping  in  a 
One-Dimensional  Aeroelastic  Problem,  Computational  Fluid  Dynamics,  6,  pp.  157-174,  1996. 


9 


Publications  Related  to  This  Research 


C.  Hindman,  M.  Balas,  and  M.  Lesoinne, “Exponential  Stability  of  Controllers  for  Fluid- 
Structure  Interaction  Using  Reduced  Order  System  Models  with  Residual  Mode  Filters”, 
Proceedings  of  OPTI  2005,  Skiathos,  Greece,  June  2005.  This  paper  uses  the  ROM/RMF 
structure  on  an  actuated  piston.  We  prove  that  exponential  stability  can  be  restored  using 
RMFs  to  augment  the  ROM-based  controller  with  actuator  dynamics. 


M.  Balas,  C.  Hindman  and  P.  Waters, “Control-Based  Reduced  Order  Structures  Models 
for  Use  in  AeroServoElastic  Simulations”,  AIAA  DTE  Summit,  Woodland  Hills,  CA, 
Nov  2004. This  summarizes  our  research  on  the  ROM/RMF  approach  but  without 
actuator  dynamics. 


Hindman,  C.,  Balas,  M.  and  Lesoinne,  M.,  Exponential  Stability  of  Controllers  Using 
Reduced  Order  Actuator  and  System  Models  with  Residual  Mode  Filters,  2000  American 
Control  Conference,  Chicago,  IL,  2000.This  develops  a  ROM/RMF  approach  for  aero- 
servoelastic  problems  but  without  stability  proofs.  The  NACA  3-d  wing  is  used  as  an 
application  here. 

Lesoinne,  M.  and  Hindman,  C.,  A  Mixed  Model  Reduction  Technique  for 
Aeroservoelasticity,  AIAA  Paper  #  99-1266,  40th  AIAA/ASCE/AHS/ASC  Structures, 
Structural  Dynamics,  and  Materials  Conference,  St.  Louis,  MO.,  Apr.  12-15,  1999.  Some 
of  the  ideas  in  this  work  were  applied  to  our  current  work  in  the  aero-elastic  problems 
above. 


10 


